Quantum Monte Carlo calculations of neutron-alpha scattering 
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We describe a new method to treat low-energy scattering problems in few-nucleon systems, and 
we apply it to the five-body case of neutron-alpha scattering. The method allows precise calculations 
of low-lying resonances and their widths. We find that a good three-nucleon interaction is crucial 
to obtain an accurate description of neutron-alpha scattering. 
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There has been significant progress recently in un- 
derstanding the ground and low-lying excited states of 
light nuclei through microscopic calculations with realis- 
tic two- and thrcc-nuclcon interactions 

dEH. These 

studies have highlighted the importance of including a 
three-nucleon interaction to obtain correct overall bind- 
ing of the p-shell nuclei p}, the ordering of states in 
10,11, 12g jgjj an( j tne stability f neutron-rich nuclei 

Modern calculations have been very successful in repro- 
ducing a large number of the experimentally-observed nu- 
clear levels up to mass 12. 

A wealth of additional experimental information is 
available in the form of low-energy scattering and reac- 
tion data. Very narrow low-energy resonances have been 
treated in the calculations discussed above as essentially 
bound states, but non-resonant scattering and broad res- 
onances have not been adequately addressed. These cases 
require the development of methods to compute scatter- 
ing states, which in turn will form the foundation for 
a quantitative, microscopic theory of low-energy nuclear 
reactions on light nuclei. Such a theory would be useful 
to astrophysics, because experimental determinations of 
crucial reaction rates are often difficult and sometimes 
impossible. 

Here we review the method previously used for quan- 
tum Monte Carlo (QMC) calculations of low-energy re- 
action properties, present an improved method, and use 
it to study low-energy neutron-alpha scattering as a five- 
body problem. We evaluate the two low- lying p-wave res- 
onances with J T =3/2 _ and l/2~, respectively, as well as 
low-energy s-wave (l/2 + ) scattering. 

In each case, we calculate the phase shift as a func- 
tion of energy for the Argonne v%s (AV18) two-nucleon 
potential [I| alone, and with either the Urbana IX (UIX) 
model Q or Illinois-2 (IL2) model [H three-nucleon po- 
tential added. We find significant differences in the calcu- 
lated phase shifts among these Hamiltonians. The AV18 
and AV18+UIX models produce too little splitting of the 
3/2~ and l/2~ states, while AV18+IL2 reproduces their 
energies and widths very well. All three models provide 
good matches to the low-energy s-wave cross sections. 



The results demonstrate that experimental phase shifts 
can be reproduced using realistic interactions with sim- 
ple modifications to the computational method used for 
bound states. They also suggest that scattering calcula- 
tions could provide sensitive constraints on models of the 
three-nucleon interaction. 

In general it is very difficult to treat quantum scat- 
tering problems with more than a very few constituents. 
Often there are many initial or final states, and it is not 
yet possible to discretize the Schrodinger equation and 
solve directly for the scattering states. Correct treatment 
of the boundary conditions for more than two outgoing 
constituents is also complicated. 

For low-energy scattering, though, there are often only 
a few two-cluster channels open. For one open channel of 
total angular momentum J and orbital angular momen- 
tum L, the wave function at large separation r of clusters 
with internal wave functions $ c i and <f> C 2 will behave as 
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where 5jl is the phase shift, k = y2/iE /h 2 , E is the 
scattering energy in the cm. frame, /i is the reduced 
mass, Fl and Gl are respectively regular and irregular 
real solutions of the Schrodinger equation with zero nu- 
clear potential, Yr, are spherical harmonics, and {• • -}j 
denotes angular momentum coupling. 

Previous quantum Monte Carlo calculations of scat- 
tering states converted the scattering problem to an 
eigenvalue problem by imposing a boundary condition 
l F(r = i?o) = at a surface radius Rq beyond the range 
of the nuclear interaction, and then solving for the en- 
ergy within this nodal surface 0,1- Eq. QJ then gives 
tan<5j£ = — Fi,(kRo) / G L(kRo) . Finding E as a function 
of i?o is then equivalent to determining Sjl(E). 

The zero boundary condition requires different Rq for 
different scattering energies. Energies near threshold re- 
quire very large Rq, which can cause numerical difficul- 
ties. For these reasons, it is preferable instead to impose 
a logarithmic derivative 7 on the wave function along the 



2 



direction n normal to the r — Rq surface: 
n • Vr* — 7* , at r = R . 



(2) 



Equation JT]) then gives 5jl from E, Rq, and 7. 

We compute nuclear energies by using two QMC 
methods as successive approximations. The variational 
Monte Carlo (VMC) method uses a trial wave func- 
tion ^t{J v ]T) specified by variational parameters that 
are adjusted to minimize the energy expectation value; 
Monte Carlo techniques are used to evaluate the integrals 
over particle position. The method and the form of ^>t 
are described in detail in Refs. [HQ. VMC calculations 
of n-a scattering were first made in the 1980s [l(j, by 
using a that explicitly went to zero at Rq. In the 
present work, correlations inside are constrained so 
that it goes to the form Eq. (JTJ at large r and satisfies 
Eq. ©. 

Green's function Monte Carlo (GFMC) takes the VMC 
trial state and evolves it in imaginary time r, 



*(r) = exp[-(i? - E )t]$>t 



(3) 



so that \&(r — ► 00) oc ^0, the lowest state with the speci- 
fied quantum numbers. The propagation time is divided 
into many small steps of length At, and Eq. ([3]) is com- 
puted by iterating the small-time-step Green's function, 



G(R',R) = (R'\exp[-(H - E )Ar]\IL) 



(4) 



acting on samples of ^t- Here R = (ri, T2, . . . , r^) is the 
spatial configuration of the nucleons, and we suppress 
spin-isospin labels for simplicity. The nuclear GFMC 
method is described in Refs. [H, Q. GFMC calculations 
with the nodal condition at Rq have been performed for 
n-a scattering previously [r"J ] and also recently for some 
atomic and condensed-matter problems 12|, |l3j . 

In this work, we restrict the GFMC simulation to clus- 
ter separations less than Ro by rejecting proposed Monte 
Carlo steps that go outside the boundary. To enforce 
the boundary condition in Eq. @, we regard the wave 
function in the restricted volume as the central region of 
a scattering wave function that fills all space, with some 
energy specified implicitly by the boundary condition. 
Then the contributions to \&(r) at the (n + l)th step in 
r come from inside and outside the r < Rq region: 

*„+i(R')= / dR cl dR c2 dr G(R',R)* n (R) 

J\r\<Rn 



I, 



|r e |>flo 



dR cl dR c2 dr e G(R, R)^„(R) . 



(5) 



Here R is written in terms of both internal cluster coordi- 
nates R c i and R C 2, and the cluster-cluster separations r 
and r e , which respectively are contained inside, or extend 
outside, of the boundary Rq. 

The contribution of the outer region can be mapped 
to an integral over the interior region by a change of 



variables, r = (Ro/r e ) 2 r e with r e — r e |, in the second 
term of Eq. ([5]). A wave function sample at the point 
R' is then the sum of one contribution from the previous 
point R and one from an "image point" at R e , located 
outside the r = Rq boundary: 

*„+i(R') = / dR c i dR c2 dr G(R', R) 

J\r\<Ro 



'\r\<Ro 



(6) 



Because G(R', R) has short range, we can use the bound- 
ary condition to obtain \t n (R e ) by linear extrapolation, 



* n (R e ) « [1 + 7 (Re - R) • &] *«(R) 



(7) 



or by fixing the energy and extrapolating with Eq. JT]). 
The GFMC simulation considers each partition of the 
nucleons into clusters with separation r. If the resulting 
|r — r e | < 1 fm, then the image contribution is used; 
otherwise the short range of G makes it insignificant. 

This method has several advantages over the nodal 
boundary condition. The most important is that Rq can 
be set to a constant value, independent of the scattering 
energy. It is also possible to calculate correlated energy 
differences in the GFMC method for different values of 
the boundary condition if they are all computed for the 
same Ro. In such cases the path integrals differ only at 
the surface, and the energy differences are more accurate 
than the individual energies. 

Neutron-alpha scattering provides a convenient test of 
several important properties of low-energy nuclear scat- 
tering. There is no bound state in the A— 5 system, but 
the 3/2 ~ channel has a sharp low-energy resonance. The 
1/2 - state is broader and higher in energy, and the com- 
bination of these two states provides a simple measure of 
spin-orbit splitting. Since the alpha particle is so tightly 
bound, simple single-channel scattering continues up to 
fairly high energies. 

The quantity most naturally given by a GFMC calcu- 
lation is the total energy of a system, and a precision in 
the vicinity of 1% has been the goal of our past calcu- 
lations. However, the quantity used to compute phase 
shifts is the energy relative to threshold, so a precision 
of 100 keV in this quantity is ^0.3% of the total energy 
in the 5 He problem. To attain this, we must pay close 
attention to the choice of Rq, construct a VPy as close as 
possible to the desired GFMC solution, and use a less 
stringent GFMC path constraint. 

First, the only a priori constraint on Rq is that it 
should be "beyond the range" of the nuclear interaction. 
We find that our GFMC result depends on Ro at the 

level of about 100 keV (out of a total of 28 MeV) 

in going from Rq = 7 fm to 9 fm. We find no further 
change going from 9 fm to 10 fm; the highest attainable 
energy (corresponding to the state with a node at the 
surface) decreases as Rq increases, so we choose Rq = 9 
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FIG. 1: (Color online) Phase shifts for n-a scattering. Filled 
symbols (with statistical errors smaller than the symbols) are 
GFMC results; dashed curves are fits described in the text; 
and solid curves are from an 7?-matrix fit to data [13 ]. 



fm. Future calculations extending to energies beyond 
this maximum-energy state should be analogous to pre- 
vious calculations of multiple bound states with the same 
quantum numbers [l5j |. 

Second, the GFMC energy also depends somewhat on 
the input We find it important to adjust pair correla- 
tions between particles in different clusters (between the 
n and constituents of the a in this case) so that Eq. ((T|) 
is enforced at large cluster separation [l6[. We also ad- 
just a parameter in that corresponds to the scattering 
energy until it matches the final GFMC energy; this typ- 
ically takes one or two iterations of the VMC and GFMC 
calculations to obtain a self-consistent result. 

Finally, in all of our A > 4 GFMC calculations, we 
use a path constraint [l| on the GFMC walk to mitigate 
the Fermion sign problem; we compute energy samples 
only after releasing the constraint for some number of 
steps to avoid biasing the results. We find that stable 
results in our scattering calculations require the use of 
80 unconstrained steps rather than the usual 20 to 40. 

In Fig. Q] we present phase shifts for all channels, com- 
puted with three different interaction models. In each 
case the AV18 potential is used as the two-nucleon inter- 
action; in the second (third) case the UIX (IL2) three- 
nucleon potential is added. We also show partial-wave 
total cross sections for the AV18+IL2 case in Fig. [3J We 
compare these results with those from a multi-channel R- 
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FIG. 2: (Color online) Calculated and ii-matrix partial-wave 
cross sections. The calculations, shown with their Monte 
Carlo error bars, are for the AV18+IL2 Hamiltonian. Stars 
show the pole energies in 3/2~ scattering for the 7?-matrix 
fit and for AV18+IL2, with the bars indicating the imaginary 
part. 



matrix analysis of the 5 He system 14| that characterizes 
the measured scattering data very well (x 2 /d.o.f. is 1.6). 
Some of the resonance parameters from that analysis are 
given in Refs. [ItJ and 18j]. Because there are more than 
2600 data points in the analysis, the uncertainties in the 
i?-matrix phase shifts are likely to be much smaller than 
the errors in the GFMC calculations. 

We have made rational polynomial fits to 
tan<5ji/fc 2i+1 , converted these to rational polyno- 
mials for the S'-matrix, and used these to find the poles 
of S. These fits are shown as dashed curves in the 
figures. For each of the two p-wave states, we find just 
one pole that is stable as the degrees of the polynomials 
are changed; we identify these as the resonance poles. 
For 3/2" the poles are at 1.19 - 0.77i, 1.39 - 0.75i, and 
0.83 - 0.35i for AV18 alone, AV18+UIX, and AV18+IL2, 
respectively, compared with 0.798 — 0.324i MeV from 
analysis of the data [l8|. The corresponding l/2~ values 
are 1.7 - 2.2i, 2.4 - 2M, and 2.3 - 2.6i MeV, compared 
with 2.07- 2. 79i MeV. The 1/2+ fits yield no stable pole, 
in agreement with the lack of a resonance in this channel 
and with the i?-matrix analysis. We have quoted all pole 
locations, and the scattering length below, so that there 
is an error of not more than 3 in the last decimal place. 

It is well known that realistic two-nucleon interactions 
alone provide insufficient spin-orbit splitting in light nu- 
clei [l9|, H3] • The figures and pole positions above con- 
firm this: although there is some spin-orbit splitting with 
AV18, it is less than half of what is needed to explain 
the data. The UIX three-nucleon potential includes a 
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two-pion-exchange term and a phenomenological short- 
range repulsion, fitted simultaneously to the triton bind- 
ing energy and the saturation density of nuclear matter. 
Adding UIX to AV18 increases the spin-orbit splitting, 
but not enough to match the data. 

In addition to the terms in the UIX model, the IL2 
potential includes a three-pion-ring exchange term and 
an additional two-pion-range term. The strengths of all 
four terms were adjusted to fit seventeen energy levels in 
nuclei up to A=8 [4j. Adding IL2 to AV18 induces even 
greater spin-orbit splitting in the fitted nuclei; the fig- 
ures and pole positions show that the experimental 5 He 
splitting is almost exactly reproduced by AV18+IL2. For 
comparison, excitation energies for this interaction have 
an RMS deviation of 600 keV from experiment for the 
states to which it was fitted. The widths of the reso- 
nances are also well reproduced by AV18+IL2. 

All three potentials produced essentially identical re- 
sults for the l/2 + case, in good agreement with the lab- 
oratory data. The s-wave scattering is dominated by a 
node in the n-a correlation at a position that roughly 
fixes the zero-energy scattering length, and the wave 
functions for all potentials contain similar structure. All 
three are consistent with a scattering length of 2.4 fm, 
compared with an experimental scattering length of 2.46 
fm Q. 

In summary, we have introduced a new QMC tech- 
nique to calculate low-energy scattering and applied it to 
neutron-alpha scattering. Models of the three- nucleon in- 
teraction that have been successful in describing bound 
levels of light nuclei also provide good descriptions of 
these scattering states. In the case of the most success- 
ful potential, AV18+IL2, the description is particularly 
good. The sensitivity of the results in Fig.[T]to the Hamil- 
tonian suggests that scattering calculations can provide 
important additional constraints on the three-nucleon in- 
teraction. We note that the energies obtained here for 
the 3/2~ pole are close to those obtained by treating it 
as a bound state while the pseudo-bound-state cal- 
culations of the broad l/2~ state proved less reliable. 

Many extensions and applications of the computa- 
tional method presented here are possible. Because 
these calculations are done in coordinate space, includ- 
ing Coulomb interactions between the clusters will pose 
no problems; the Ft and Gt of Eq. ([1]) are then simply 
Coulomb, rather than Bessel, functions. Extensions to 
nucleus-nucleus (rather than nucleon-nucleus) scattering 
should also be feasible. Besides pure scattering problems 
and the use of scattering wave functions to compute elec- 
troweak cross sections [21| and hadronic parity violation, 
the scattering method should be applicable to cases with 
more than one open two-cluster channel. For A c h such 
channels, a set of jV c h linearly independent solutions of 
the Schrodinger equation at the same energy determines 
the iS-matrix. The solutions could be obtained either by 
varying the boundary conditions at the surface, or by 



calculating the derivative of the energy with respect to 
changes in the logarithmic derivatives. 
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